Transverse Radial Expansion and Directed Flow 

Sergei A. Voloshin* 

Physikalisches Institut der Universitdt Heidelberg, Heidelberg, Germany 

(February 9, 2008) 
The effects of an interplay of radial expansion of the thermalized system 

created in a heavy ion collision and directed flow are discussed. It is shown 

v^Q , that the study of azimuthal anisotropy of particle distribution as a function 

On 

ON ' of rapidity and transverse momentum could reveal important information on 



> 

o 



both radial and directed flow. 

Recently, the study of collective flow in nuclear collisions at high energies has attracted 
an increased attention of both theoreticians and experimentalists. There are several reasons 
for that: i) the observation of anisotropic flow at the AGS [|T],3 and, probably, at the 
SPS II energies, ii) better theoretical understanding of the relation between appearance 



and development of flow pattern during the fireball evolution and the processes such as 
thermalization, creation of quark-gluon plasma, phase transitions, etc. ^ ne study 

of mean field effects [|]J§, iv) the importance of flow for other measurements such as two- 



particle interferometry f9Hl"3|, v) development of new techniques suitable for flow study at 
high energies p^JT5[] . Although all forms of flow are interrelated and represent only different 
parts of one global picture, usually people discuss different forms of collective flow, such as 
longitudinal expansion, radial transverse expansion, directed flow, elliptic flow |15|||. 

Flow introduces strong space-momentum correlation in the particle production. Particles 
with a given rapidity and transverse momentum are produced only by some part of the entire 
source. This part we call as an effective source. It moves with the rapidity close to that 
of the particles [§,[[]]] in the midrapidity region, or slightly less than that for particles in 



the fragmentation region [ 13| . The subject of the current study is transverse (radial and 



directed) flow, namely the questions how azimuthally symmetric radial flow interferes with 
directed flow and affects the directed flow signal. 
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It is convenient to characterize directed flow by v% , the amplitude of the first harmonic in 
the Fourier decomposition of the particle azimuthal distribution [|14|| . Depending on which 
particular distribution is studied, the coefficients in the Fourier decomposition depend on 
different variables. If it is a 2-dimensional rapidity and azimuthal angle distribution, then v 1 
depends only on rapidity, and vi = (p x ) / (pt), where (pt) is the mean transverse momentum, 
and (p x ) is an another widely used quantity to describe flow, the mean projection of the 
transverse momentum onto the reaction plane. If one studies the production of particles 
with a given rapidity and transverse momentum, then 

d 3 N d 3 N 1 d 2 N , . , A . (A 

£ W = = + 2t " c ° sM + 2vM ) + •■■)• (1) 

and v\ depends on two variables y and pt- The dependence t>i(pt) for particles with a given 
rapidity is the main topic of the current study. We will show that the radial expansion of 
the source results in specific shapes of this function, which can be studied experimentally, 
providing an information on the directed as well as radial components of the effective source 
velocity. 



First studies of event shape anisotropy as a function of transverse momentum flqjr7| have 
shown very interesting results. Such a study using Fourier decomposition of the azimuthal 
distribution promises better quantitative description of the effect. One of the advantages 
would be the possibility to correct the results for the reaction plane resolution [Hjni]. The 
analysis of v i{p t ) is to a large extend independent from the uncertainties in the p t dependence 
of spectrometer acceptance and efficiency, and in this sense such an analysis has preference 
in comparison with the analysis of triple differential distributions. 

Transverse directed flow is a result of a movement of an effective source in the transverse 
plane (below we assume that this movement is along the "x" axis). In the source rest frame 
the first moment of azimuthal distribution (vi) is zero. The final anisotropy (here and 
below we discuss only the first harmonic of the azimuthal distribution) appears only as a 
consequence of the source movement in the transverse direction. Below we derive a general 
expression for v±, assuming that the invariant distribution in the source rest frame is known 
and (for simplicity) azimuthally symmetric: 

d 3 N 

E— = J(p t ,y) = J(p t ,y). (2) 



Source movement in the transverse direction results in J — > J' = J(jp t ',y), where pt' is 
defined by the Lorentz shift with velocity (3 a along the "x" axis. Assuming that the shift is 
small (the directed flow velocity expected to be f3 a < 0.1), p' x = p x — (3 a E, and taking into 
account that dp t /dp x = cos (</>), one gets: 

= J d<j)J(pt',y) cos(0) _ f3 a E din J (p t , y) 
Vl fd<f>J(p t >,y) 2 dpt ' U 

A comparison of exact numerical calculations for the models discussed below with the results 

of application of this formula shows that for the values of (3 a < 0.1, the formula (§) is accurate 

at the level of a few percent. 

Below we apply the derived formula to a few particular cases, when the thermal source 
undergoes an isotropic expansion. We first consider an expansion in a non-relativistic case, 
then go to a relativistic generalization, and, finally, discuss a model of isotropically expanding 
thermal shell |18|], widely used in the analysis of BEVALAC and SIS data |l9|p0| . Using 
simple models permits us to perform all calculations analytically to keep clear the qualitative 
features of the effect. Note that directed as well as isotropic expansion velocities are defined 
in the frame moving longitudinally with the rapidity of the effective source. 

Directed flow of protons is the most pronounced and less affected by another effects 
such as shadowing. Protons are relatively heavy particles ( m > T, m is the proton 
mass, T is the temperature), and often can be treated as non-relativistic particles. This is 
certainly justified in the source rest frame, where the proton kinetic energy is of the order 
of temperature E* — m ~ T <C m. To have possibility to treat protons non-relativistically 
in the transverse direction in the analysis frame we have to assume that m t — m<^m, what 
restricts the non-relativistic consideration to the region of p% < 0.5 GeV. 

The transverse isotropic expansion of the source can be described as a superposition of 
different sources moving radially with an expansion velocity (3q. Then the (non-relativistic) 
transverse momentum distribution of protons from a radially expanding thermal source can 
be written as: 

1 d 2 N _ 1 fdij>exp ( - P° cos W) 2 + (Pv Z Po sin Wn / 4 n 



N dp t (2tt) 2 (2mT) J T r \ 2mT 
where po = m/5 . The integration over tp (the orientation of the expansion velocity) results 
in the distribution: 
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where I is the modified Bessel function, and £ = j3 p t /T. 
Using the formula (|3]) one gets an expression for v\\ 

In our analysis we work in the frame moving longitudinally with the effective source rapidity. 
In this frame the longitudinal particle momenta usually can be neglected, what was done in 
the derivation above. If one considers the particle production with rapidities far from the 
rapidity of the effective source, when particle longitudinal momenta are large, one should 



make a substitution m — > ym 2 + pf. 

It follows from ([]) that in the case without radial expansion v\ linearly depends on 
the transverse momentum with the slope of 0.5/3 a /T. The radial expansion of the system 
decreases the directed flow signal changing also the shape of Vi(p t ) in the low p t region. 
Remarkably, that at some parameter values, namely such that m(5\j2 > T (for T = 100 MeV 
it corresponds to (3q > 0.46), Vi(p t ) has a region (at small p t ), where it is negative (see Fig. 1). 
Physically it corresponds to the case, when particle production with such a value of p t is 
more probable from the part of the effective source which moves to the opposite direction 
than the flow direction (in this case directed flow and "expansion" flow compensate each 
other). 

The relativistic generalization of the above formulae is strait-forward, although results 
into somewhat less transparent expressions. The invariant distribution in this case can be 
written in the form: 

Joe J d^E*e~ E 'l T 

= J d^(E cosh(y t ) - p t cos(^) sinh(y t )) exp(-(E cosh(y t ) - p t cos(^) smh(y t ))/T) 

= Te-*(xJo(0-tfi(0). (7) 

where \ — E cosh(y t )/T, £ = p t sinh(y t )/T; y t = 0.51n((l + Po)/(l — A))) is the transverse 
rapidity; J and I\ are modified Bessel functions. Such a distribution yields: 

h{g)E? sinh(^) cosh(^) + J (e)(p t Tcosh(^) - ETsmh(y t ) 2 ) 



f3 a p t cosh(y t ) 
Vl = 



p t Tcosh(y t )( X I (0-Zh(0) 



■ (8) 



Note, that the energy E is measured in the frame moving longitudinally with the same 
velocity as an effective source. Studying particle production close to the fragmentation 
region, one should remember, that the rapidity of the effective source could be less than 
that of the particles In this case E = m t cosh(?/* — y), where y* is the effective source 
rapidity 

If the particles are emitted from a thermalized spherical shell at temperature T, ex- 
panding with velocity (3 r (we use different notation to distinguish (3 r from /3 , the expansion 
velocity in the transverse plane), the expected invariant distribution has the form [|TB| , |iT? 

sinh(a 



J oc exp(— E'y/T) 



a 



■(E + T/i) -cosh(a)T/7 



(9) 



where 7 = 1/y 1 — (3f, and a = (3 r p"y/T. An application of the prescription (|3|) gives: 

V\ = ^2^? sinh(a;) /a — cosh(a)T/ 7 + \{E + T j 7) (cosh(a) /a — sinh(a;) / a 2 ) 

- smh(a)T/j)p r E/p} / {(E + T/j) sinh(a)/a - cosh(a)T/ 7 } (10) 

Qualitatively this model predicts the same dependencies Ui(p*) as the model with the ex- 
pansion only in the transverse plane. But the value of the expansion velocity needed in this 
model for v\{p t ) exhibits a dip at low p t is about 15% larger due to the 3-dimensional type 
of the expansion. 

Summarizing, it is shown that a study of the dependence of the amplitude of the first 
harmonic in the azimuthal distributions on transverse momentum reveals important features 
of directed and transverse radial flow. The magnitude and the shape of Vi(p t ) is very 
sensitive to such parameters of the system as an effective temperature, radial expansion 



and directed flow velocities. Strong radial expansion (approximately (3q > y2T/m) would 
evidence itself in the directed flow signal showing an anti-flow of the particles with low 
transverse momentum. The preliminary results of E877 Collaboration presented at HIP- 
AGS '96 Workshop f2]J suggests that such a regime could be the case in gold-gold collisions 
at the AGS. 

The author acknowledge the discussions with W.-C. Chang. 
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I. FIGURE CAPTION 

FIG. 1. vxipt): solid lines - Eq. (§), dashed lines - Eq. (§). T = 0.1 GeV, (3 a 
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FIG. 1. «i(p t ): solid lines - Eq. (6), dashed lines - Eq. (8). T = 0.1 GeV, f3 a = 0.1. 
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